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Abstract 

In the Landau-de Gennes theoretical framework of a Q-tensor descrip- 
tion of nematic liquid crystals, we consider a radial hedgehog defect with 
strong anchoring conditions in a ball B C I s . We show that the scalar 
order parameter is monotonic, and we prove uniqueness of the minimizing 
hedgehog below the spinodal temperature T*. 

1 Introduction and notations 

In the theory of Landau and de Gennes, a nematic liquid crystal is described 
by its tensor order parameter, the so-called Q-tensor [2] j which is a traceless 
symmetric matrix : denoting by S3 (R) the set of symmetric 3x3 matrices with 
real entries, it holds 

QeS 3 (R), trQ = 0. 

A null Q-tensor corresponds to an isotropic liquid, a Q-tensor with two equal 
(non zero) eigenvalues correspond to a uniaxial nematic, and a Q-tensor with 
three distinct eigenvalues corresponds to a biaxial nematic [2]. 

In the case of a uniaxial nematic, the Q-tensor can be written 

Q = s(n (gin- ~1), (1) 
o 

where n is a unit vector and s£K. The unit vector n is called the director of the 
nematic : it corresponds to the preferred direction, along which the molecules 
tend to align. The scalar s is called the scalar order parameter : it quantifies 
the degree of alignment of the molecules along this director Jj . The value s = 
corresponds to the isotropic state. 
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We are going to consider a nematic crystal liquid filling a cavity f2 C R 3 . The 
system is spatially inhomogeneous : the tensor order parameter Q(x) depends 
on x € and at equilibrium it should minimize the energy 

E(Q)= [ MQ) + fe(Q) dx, 

where 

MQ) = iatr(Q 2 ) - htr(Q 3 ) + \cMQ 2 )]\ (2) 

and 

/e(Q)-§|VQ| 2 . 

The bulk energy density fb describes the nematic-isotropic transition in a spa- 
tially homogeneous system, while f e accounts for elastic deformations of the 
nematic. Here we work with a particularly simple form of the elastic energy 
density / e , using the so-called one-constant approximation [2]. For the bulk 
energy density fb, terms of higher order could be taken into account, but the 
expansion up to fourth order is usually considered - because it is the simplest 
one that allows to describe the transition 2 . The constants a, b, c and L are 
material- and temperature-dependent, and in practice only a depends signifi- 
cantly on the temperature. This dependence is usually assumed to be affine 

a = a (T - T*), 

where T* is the so-called spinodal temperature. 

It can be shown [51 Proposition 1] that, for spatially homogeneous systems, 
the Q tensors minimizing the bulk energy density fb are necessarily uniaxial (or 
isotropic), i.e. of the form ([T]). 

For uniaxial Q-tensors of the form (JTJ, fb depends only on the scalar order 
parameter s : 

MQ) = Ms) = ^ - + 

One can distinguish between different temperature regimes, with three crit- 
ical temperatures T* <T C < T+. 

• In the temperature regime where T < T*, i.e. a < 0, the isotropic state is 
unstable : fb(s) has no local minimum at s — 0, and has a global minimum 
at some value s + > 0. 

• For T* < T < T c , M s ) s ^ nas a global minimum at some s + > 0, but 
also has a local minimum at s = 0. 

• For T c < T < T + , fb(s) still has a local minimum at some s+ > 0, but 
the global minimum is at s = 0. 

• Above T + the only critical point of fb(s) is at s = 0. 
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We always consider the regime where there is a positive global minimizer of 
fb{s) (i.e. T < T c ) : the nematic state is stable (but the isotropic state also). 
For the uniqueness result (Theorem 0| we will further restrict ourselves to the 
case where there is no local minimizer at s — (i.e. T < T*) : the isotropic 
state is not stable. 

Here we study a nematic droplet : the container is 

SI = B R = |i£l 3 : |i| < Rj . 

A proper rescaling 3 allows to work with dimensionless variables and to get 
rid of the material-dependent constants : we can take the energy to be 

e{q)= L (k\ vQ \ 2+h{Q) ) dx ' (3) 

with the rescaled bulk energy density 

fb(Q) = \ MQ 2 ) - V6 tr(Q 3 ) + 1 tr(Q 2 ) 2 + C(t). 

Here t — 27ac/b 2 = 27ao(T — T*)c/b 2 is a dimensionless temperature and 
C(t) is a constant chosen so that min/b = 0. 

This non-dimensionalization of the problem is described in [X] below. 

In terms of the reduced temperature t, the critical temperatures are < 1 < 
9/8 [3], and the different temperature regimes mentioned above are: 

• For t < 0, the isotropic state is not a local minimizer of /(,. (This corre- 
sponds to T < T*.) 

• For < t < 1, the isotropic state is a local minimizer of /& but the nematic 
state is the global minimizer. (This corresponds to T* < T < T c .) 

• For 1 < t < 9/8 the global minimizer of is the isotropic state. However, 
the nematic state is a local minimizer. (This corresponds to T c < T < T + .) 

• For t > 9/8 the nematic state is not defined anymore, since the only local 
mimimizer is isotropic. (This corresponds to T > T + .) 

Hence in the following we assume t < 1, and for the uniqueness result (Theo- 
rem|5) t < 0. 

We restrict ourselves to the case of strong radial anchoring, which means 
that we have the Dirichlet condition 

Q^ = \^Mw\^¥\--^ forN = i? ' (4) 

where /i+ is the scalar order parameter minimizing fb(Q) — if and only if 
Q is uniaxial with scalar order parameter ^3/2h+ . Its value is given by [3] 

3 + V9^8l 
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The admissible space for Q is then 

W Q = {QE W 1 ' 2 (B R ,S 3 {R)),trQ = 0, Q satisfies gj)} . 

In general the minimizing Q-tensor may not exhibit complete spherical sym- 
metry: symmetry breaking configurations have been observed experimentally 
and theoretically (see [5] and references therein). However, here we restrict 
ourselves to the study of spherically symmetric critical points of E. 

To this end we define an action of the special orthogonal group SO (3) on 
the space of admissible Q-tensor configurations Wq, in the following way: 

(U -Q){x) = UQ( t Ux) t U, U G SO(3). 

The spherically symmetric configurations, i.e. those Q G Wq satisfying 

U ■Q = Q VU G SO(3), 

are precisely the maps of the form 

0^) = ^(W)(r®R-5 1 )' (5) 

where the scalar order parameter h belongs to the admissible space 

W h = {he L 2 (0,R;dr),ti G L 2 (0,R;r 2 dr),h{R) = h+) . 

An equilibrium configuration described by such Q is called radial hedgehog. 
Hence we consider, for h G Wh, the energy functional 




where 

9(h) = MQh) = \h 2 -h 3 + h 4 + C(t). 

The main goal of this article is to study minimizers of this functional, and in 
particular their monotonicity and uniqueness. 

The problem of studying droplets of nematic liquid crystals is of great phys- 
ical relevance. It plays an important role in some electro-optic applications, 
like polymer dispersed liquid crystals (PDLC) devices (see the review article [7] 
and references therein). Moreover, defects in general, and the radial hedgehog 
defect in particular, usually exhibit universal features [5], which makes it im- 
portant to understand their structure better. The monotonicity of the scalar 
order parameter, which is proved in this article, is physically intuitive and had 
been observed numerically [TT], but no general mathematical proof of it was 
available. The uniqueness of this minimizing profile constitutes the main result 
of the present article. 

The spherically symmetric critical points of E have already been studied in 
a recent paper by A. Majumdar [9]. A. Majumdar also addresses the questions 
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of uniqueness and monotonicity of the scalar order parameter h, in the limiting 
cases of low temperature (t — > — oo) and big radius (R — > oo), using techniques 
based on analogies with Ginzburg-Landau vortices. In the present paper these 
questions are studied without any great restriction on temperature and radius, 
using elementary properties of the energy functional I, and the differential equa- 
tion its critical points satisfy. 

It should be mentioned here that there exist other continuum theories for 
nematic liquid crystals. The two most important are the Oseen-Frank model and 
the Ericksen model, which can both be viewed as simplifications of the Landau- 
de Gennes theory. In particular they only account for uniaxial configurations. 
For more information on the mathematical theory of liquid crystal, see the 
review article [6] and the references therein. 

The paper is organized as follows. In Section[2j we recall a few basic proper- 
ties of minimizers of I. Among them is a singular ordinary differential equation 
(0 satisfied by the scalar order parameter. In Section [3] we will see that, even 
though this equation ([7J is singular, it leads to a well-defined Cauchy problem. 
In Section 2] we show that h is monotonic. And in Section [5] we prove our main 
result, Theorem [4] asserting that I has a unique minimizer (in the temperature 
regime T <T*). 

2 First properties 

The following proposition summarizes the first basic properties of the radial 
hedgehog. It is very similar to [HI Proposition 2.1.], but we give a proof here for 
completeness. It strongly relies on the symmetries of the problem. 

Proposition 1. (i) There exists a minimizer h € Wh of 1(h). It satisfies the 
singular ordinary differential equation 

h"(r) + -h\r)-^P-=g\h(r)), re(0,R). (7) 
r r z 

(ii) The corresponding Q-tensor Qh € Wq (as defined by above) is a crit- 
ical point of the non-restricted energy functional E, h is analytical on 
[0, R), and it holds 

h(0) = h'(0) = 0. 

(Hi) For every r G [0, R], 

< h(r) < h+. 

Proof. Let (h n ) be a minimizing sequence for I : 

h n E W h , I(h n ) — > mini. 

Define Q n = Qh„ € Wq. Since 

/ |VQ„| 2 dx < 2E(Q n ) = &TTl(h n ), 
JB R 
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and Q n has fixed boundary conditions independent of n, the sequence (Q n ) 
is bounded in the Hilbert space W := W 1 ' 2 (Br ) S 3 (R)), so, up to taking a 
subsequence, it converges weakly to some Q £ W. 

Since the integrand defining the energy functional E is convex in VQ, E is 
weakly lower semi-continuous (see e.g. [T| Theorem 1.3]), so that 

E(Q) < liminf E(Q n ) = 47rliminf I(h n ) = 4ttMI. 

The set of spherically symmetric Q-tensors is weakly closed in W because it 
is a closed subspace of W. Therefore the weak limit Q has to be a spherically 
symmetric Q-tensor: there exists h £ Wh such that Q = Q/j. It holds 

/(/,) = J_ S (Q) < inf /, 

hence h is a minimizer of 1(h). 

Classically, the fact that h is a critical point of / provides us with equation 
([7]): considering, for small t a perturbation h + t(p, where ip £ C£°(0,R), the 
vanishing of the first variation reads 




r 2 (h'ip 1 + —hip + g'(h)ip) dr = 0, 



which becomes, after an integration by part, 




( - r 2 h" - 2rti + 6r 2 h + r 2 g{h))<pdr = 0, 



therefore (0 holds in the sense of distributions in (0,R), and hence it holds 
almost everywhere. This proves (i). 

To prove (n), we use the Principle of Symmetric Criticality |13j . For clarity, 
we state this principle in the form we will use here : if a group G acts linearly 
and isometrically on a Hilbert space H, and E is the set of points fixed by this 
action, then for any G-invariant functional / £ C 1 (i?, K), a critical point ieS 
of /|£ is a critical point of the non-restricted functional /. 

The action of SO(3) defined above is linear and isometric. Our energy 
functional E is invariant under this action, and it is C 1 because fb is C , 
|/ b (0)|<c|Q| 6 and|^(Q)|<c|Q| 5 . 

Hence the Principle of Symmetric Criticality applies, and Qh is a critical 
point of the energy functional E : it satisfies the Euler-Lagrange equations 

AQij = t Q y - 3V6 (Q lk Q k] - -Sij tr(Q 2 )) + 2 Q l3 tr(Q 2 ). (8) 

Elliptic regularity theory (T2J Theorem 6.7.6] tells us that Qh is analytical on 
-Bfl, so that h is analytical on [0, R), since it holds, for any unit vector uq £ K 3 , 

h(r) = ^u -Q(ru Q )u . (9) 
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As a consequence, it must hold 

h(0) = 0, 



since 



\VQ\ 2 (x) = h'(\x\f + 



6fo(M) 2 
\x\* 



is continuous. And it must further hold 



h'(0) = 0, 



since formula ((9|) extends h to a smooth and even function. It concludes the 
proof of (ii). 

The analyticity of h away from the origin was also clear from the ordinary 
differential equation ([7]) it satisfies (and also at r = R with this argument), but 
regularity at the origin was not obvious, since functions in the admissible space 
Wh may even not have a finite limit at the origin. 

As to (Hi), the upper bound follows from a maximum principle |101 Propo- 
sition 3] for equation © satisfied by Qh- The lower bound follows from 



since 1(h) is minimal. Hence, the set {h < 0} is empty, and h > as claimed. □ 



We have seen that solutions of ([7]) which lie in Wh must verify h(0) = h'(0) = 0. 
Their first degree of freedom is h"(0). As a matter of fact, one can rewrite 
equation (UJ) as 



and use the following lemma to show the existence of a unique solution of (J7J 
with fixed second derivative at the origin. 

Even though the ordinary differential equation ([Jj is singular at the origin, 
Lemma[5]shows that this singularity is, in some sense, just a pseudo-singularity: 
there is still a well-defined Cauchy problem at the origin, and the proof is similar 
to the non singular case (Picard's Existence Theorem). 

Lemma 2. Let a, j3 € K, a < 1. Define 7 := 1 — a — f3 and assume 7 > 0. Let 
F : K — > K be locally Lipschitz with F(0) = 0. Let a G R. There exists some 
R a > and a function h a : (0, R a ) — > K ; unique solution of 




3 The ODE satisfied by h 



d_ - 

dr L 



r 



dr 



(r 3 h(r)) =r- 1 g'(h(r)) 1 r£(0,R), 



— r — 

ir I dr 



[r^h(r)] 



r a+p F(h(r)), 



Vr G (0, R a ) 



(10) 




a 



as r — >• 0. 



(11) 
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Proof. It suffices to show existence and uniqueness in a neighbourhood of the 
origin. 

Let us first remark that the problem (I10l) - (|lip is equivalent, near zero, to 
the integral equation 

h(r) = ar 1 + r _/3 f p~ Q f t a+f) F(h{t)) dtdp Vr e (0, e), (12) 
Jo Jo 

for a continuous h verifying 

\h{r)\<Cr' 1 Vr€(0,e). (13) 

Indeed, if h satisfies the boundary condition (|TTj). then there exists e and C 
such that h verifies (fT3")) . If, furthermore, h satisfies (|10[) . then an integration of 
JTDJl yields 

^- [r^/il = Ar-" + r~ Q / t a+fs F(K) dt. 
dr I J 7 

A second integration (which is possible since a < 1) leads to: 

fcfr) = Br"' 3 + — ^— r 7 + r-* 3 / p" a f t a+p F(h) dtdp. 
1 - a Jo Jo 

The boundary condition (|11|) imposes B = (because — /3 < 7 as a < 1) and 
A = a(l - a). 

Conversely, if a continuous h verifies (TT2")) and (TT3")) . then derivating twice 
yields IjlOp. Moreover the second term in the right-hand side of (fT2]) is negligible 
compared to ar 1 so that the boundary condition (fTTj) holds. 

Let .R' be the Lipschitz constant of F on [— |o| — 1, |o| + 1]. Let 



' 2(3 -a) 
(l + \a\)K 



and e = min(<5, 1). We define a sequence of functions h n : (0,e] — > M in the 
following way : 

ft»+i(r) = ar~<+r-P f p~ a f t a+0 F(h n (t)) dtdp, (14) 



and ho(r) — ar 1 . The sequence (h n ) is well-defined if for all n, 

\h n (r)\<(l + \a\y Vre(0,e]. (15) 

Indeed, (fT5|) implies the convergence of the integral in (fT4"|) . We prove (TT5|) by 
induction. Suppose (1151) is true for some n. Then, making use of 

\F(y)\ = \F(y) ~ F(0)\ < K\y\ for \y\ < 1 + \a\ 

we obtain, using ()14[) and the induction assumption, we find that 

\h n+1 (r) -ar~<\< K(l + \a\) r - f3 £ p~ a £ ' t a+ ^f dtdp = ^ . 
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This completes the induction since r 2 /S 2 < 1 for r e (0, e). 

We next show that, for every r € (0, e] and every bGN, we have 

\K +1 {r) - h n (r)\ < (1 + \a\) ■ , (16) 

nl 

where C = K/(6- 2a). 

This follows by induction via the following calculation : 

\h n+1 (r) - h n (r)\ < r- 1 * [ p~ a T t a+f3 \F(h n (t)) - F(h n ^(t))\ dtdp 



o jo 



< Kr- p f p- a ( P t a+fi \h n (t)-h n -i(t)\dtdp 
Jo Jo 

^n-l rr rp 



< 



(1 + \a\)K^-—r-P / p- a / t a +H 2n - 2 ^ dtdp 
(n-1)! J a J 



(n - 1)! 2n(2n - a + 1) 

fm 7+2n 

< (i + H) — . 

n! 

Inequality (|16[) ensures that the sequence (h n ) converges uniformly on (0, e], 
to some continuous function h. The domination |/i„(r)| < (1 + lal)?" 7 allows to 
take the limit for n — > oo in (fl4| . so that h verifies the integral equation (fT2|) : 
ft, is a solution to problem (fTOj) - (jTT |) . 

To show uniqueness, assume h is a solution of pH)) - ([TT|) . Then for some e it 
holds 

\h(r) - ar 7 | < r 7 Vre(0,e), 
and (|12p . An induction similar to the one performed above allows to show 



\h(r)-h n (r)\ <(1 + M) r , Vre(0,e), 



so that h = h and the problem has a unique solution. □ 

Taking a = — 4, /? = 3 and thus 7 = 2, Lemma [2] directly applies to our 
equation, yielding the result claimed above. 

4 Monotonicity of h 

Here we keep assuming that we are in the temperature regime T < T c (i.e. 
t < 1). 

Proposition 3. ///i£ Wh is a local minimizer of I and satisfies < h < 
then h is increasing on (0,i?). In addition, we have 

h'>0 in (0,R]. 
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Proof. First we will show that 

ti>0 in (0,R). 

Indeed, let 

W hfi = {ip € L 2 (0, R; dr),ip' € i 2 (0, R; r 2 dr), ip(R) = 0} . 
For ijj € Wh,o> we have, as s — > 0, 

I(h + tn/>) = 1(h) + sl'(h) ■ i, + S —I"{h) • (V>, rf>) + o(s 2 ), (17) 

with r 

I'(h)-i>= f r 2 (/iV' + 6^+g'(h)ip)dr 
Jo r A 

and 

I"(h)-{il>,il>)= I r 2 (^ 2 + 6^+g"(h)^ 2 )dr. 
Jo r 

Validity of (fl7|) follows easily from the fact that g is polynomial and h is 
bounded. Since h is a local minimizer, we have I'(h) ■ ip = (this is the Euler 
Lagrange equation) and /"(/i) • > for every i/j e W/i,o- 

Our analysis of the Cauchy problem in Section [3] implies that h' > on (0, <5) 
for some positive 5. Let us assume that h is not monotonic in (0, R): there exist 
< Ri < R 2 < R such that h' = at i?i and R 2 and ft' < on (i? l5 R 2 ). 

Let t/i = h'l( RuR2 j € Wh,o- We are going to show that I"(h) ■ (ip,i/>) < 0, 
which will contradict the fact that h is a local minimizer. 

Differentiating once the equation |7} satisfied by h, one gets 

ft< 3 > + -h" - \ti + = g"(h)h'. (18) 

Multiplying (TT51) by r 2 /i' and integrating from Ri to i?2, we obtain, after an 
integration by parts (using h'(R\) — h'(R2) = 0), 

[ R2 ( - r 2 (h") 2 - 8(h') 2 + —hh') dr = r 2 g"(h)(h') 2 dr. 
Jr 1 V r ' JR! 

This can be rewritten as 

I"{h)-(ip,^)= [ 2 —hh'dr-2 [ 2 {h') 2 dr. 

We obtain the desired contradiction since h > and hi < in (Ri,R 2 ). 
Eventually we can actually show that 

h'>0 in (0,-R]. 
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Assume first that for some ro G (0, R) we have h'(ro) = then since h! > 
it must hold h"(r ) = 0. Therefore, using equation (fT8j) we get 

^ (3) M = -™h(ro) < 0, 
r o 

contradicting h' >0. 

It remains to show > 0. Assume h!(R) = 0. Using equation 

satisfied by ft, together with the fact that 

g'(h(R)) = g'(h + ) = 0, 

we have 

h"(R) = 6^ > 0, 
again contradicting h' > 0. □ 

5 Uniqueness of h 

In this section we assume t < 0, so that the function g is decreasing in [0, h+\. 

Theorem 4. In i/ie temperature regime T < T* , the energy junctional I admits 
a unique global minimizer h € Wh- 



Proof. The Euler Lagrange equations ([5]). may be rewritten as 

AQy = -J^-(Q) + % tr(Q 2 ). (19) 
Solutions of (|19p satisfy the following Pohozaev identity [TD1 Lemma 2]: 



i / |z • VQ| 2 ^(x) - | / |VQ| 2 d^ar) + \ [ |VQ| 2 dx = -<3 f j 

U JOBr 1 JdBR 1 J Br J Br 

where da stands for the usual measure on the sphere OBr. 

Using the definition of the energy E, this identity is equivalent to 



f b (Q) dx, 



E(Q) + 2 f f b (Q) dx = | / |VQ| 2 dcr(x) - 1 

JBr 1 JOBr M 



\x ■ VQf da(x). (20) 

dB R 



Now let us assume there exist hi ^ £ W/, satisfying 
/(ft!) = /(ft 2 ) = mini. 

Applying ([20"]) to Qh t for i — 1,2 and subtracting the two equations, we find 

Stt/ r 2 ( 5 (/ il )- 5 (/ l2 ))dr = 27r(^(i?) 2 -/ l ' 1 (^) 2 )- (21) 
Jo 
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The two functions hi and hi are distinct solutions of the Euler Lagrange 
equation ([7]) corresponding to the minimization of /. Hence Lemma [2] ensures 
that they cannot coincide on a neighbourhood of zero : we have for instance 
h% > hi on (0, e). Actually they can never meet before R : let us assume the 
existence of a ro € (0, R) such that hi(ro) = hi(ro). Then we obtain a new 
minimizer h of /, 

h = h 2 l(0,r o ) + hll( r0t R). 

Let us prove the fact that h is a minimizer : firstly it lies in the admissible 
space, and, denoting by e[h] the energy density (1(h) = e[h]dr), we have 

pr pR 

I(hi) < 7(/iil(o,r ) + ^2l(r ,-R)) = / e[hi] dr + / e[h 2 ]dr, 

JO Jr 

since hi is a minimizer and hilm ro \ + hii-i ro ,R) ^ es m the admissible space. 
Therefore it holds 

pro pro 

/ e[fa] dr< e[hi] dr, 
Jo Jo 

and adding e[hi]dr on both sides of this inequality exactly yields 1(h) < 
I (hi), hence proving the affirmation. 

In particular, since ft, is a minimizer it must be analytic. But at r Q all its 
right derivatives are equal to those of hi, thus h — hi on a neighbourhood of ro, 
so that hi — /12 which contradicts our primary assumption. Therefore it holds 

hi < hi on (0, R). 

In particular it implies, together with hi(R) = hi(R), that 

h[(R) > h' 2 (R), 

so the right-hand side of (f2"Tj) is non positive. 

On the other hand, since the energy density g is decreasing on [0,/i+], and 
hi < hi, the left-hand side of (|21[) is positive : we obtain a contradiction. □ 

A different but related question is whether limiting maps obtained in the 
limit R — > oo satisfy the ordinary differential equation with the limiting 
boundary condition h(oo) = h-\-. Uniqueness of solutions to this boundary 
value problem, without minimality assumptions, is established in [14J . For a 
discussion about uniqueness in the limit t — ¥ oo, see also [5J Section 3] and the 
references therein. 
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A Non-dimensionalization 

We define, as in [3] , the typical length £ and the reduced temperature t, in terms 
of which one can express a dimensionless version of the problem : 



Making the change of variable x = £x in the integral © defining the energy 
then gives 




For Q S Wq define Q by 



Q(x) = qoQ(£x) for \x\ < R 



R 



where 





Direct computation shows that 




Therefore it holds 




where the rescaled energy E is defined by 
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